FILTERING METHOD AND APPARATUS 



BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to a filtering method and 
apparatus for use in decoding the digital signal from the 
frequency domain to the time domain. 

Description of the Prior Art 

Conventionally, a method for decoding the digital signal 
intheMP3 (MPEG Audio Layer III ) or the like from the frequency 
domain to the time domain was well known that employed the 
filtering with a hybrid filter bank. 

FIG . 9 is a diagram showing the configuration of functions 
of an encoder 100 and a decoder 200 in accordance with the 
MP 3 standards. The main functions of the encoder 100 and the 
decoder 200 will be set forth below with reference to FIG. 
9. 

In FIG. 9, the encoder 100 comprises a hybrid filter bank 
(hereinafter referred to as an "HFB" ) 110, a psychoacoustic 
analysis portion 120, a iteration loop 130, a Huffman coding 
portion 140, a side information codingportion 150 , and a stream 
composing portion 160. 

The HFB 110 further comprises a sub-band analysis filter 
bank 111, an adaptive block length MDCT (Modified Discrete 
Cosine Transform) 112, and an alias reduction butterfly 113. 
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The sub-band analysis filter bank 111 cuts out an audio 
signal input into the encoder 100 by a predetermined window 
among the windows as shown in FIG. 10, and divides the audio 
signal from the time domain into a plurality of frequency 
5 domains. The adaptive block length MDCT 112 transforms each 
of these signals into more subdivided frequency spectrum in 
accordance with a predetermined transformation expression. 
The alias reduction butterfly 113 makes a predetermined 
butterfly operation on these frequency spectrum for output 

10 to the iteration loop 130. The adaptive block length MDCT 
112 makes a transformation process into the frequency spectrum 
in accordance with weighting of the data based on the analysis 
result of the psychoacoustic analysis portion 120 (i.e., a 
bit allocation based on a psychoacoustic model). 

15 The iteration loop 130 makes the non-linear quantization 

of the input frequency spectrum in accordance with weighting 
of the data based on the analysis result of the psychoacoustic 
analysis portion 120 . Also, the iteration loop 130 calculates 
a scale factor (the information regarding the weighting of 

20 the data based on the analysis result of the psychoacoustic 
analysis portion 120) from an output signal (frequency 
spectrum) oftheHFGllO, and outputs it to the side information 
coding portion 150. 

The Huffman coding portion 140 makes the Huffman coding 

25 of a signal quantized by the iteration loop 130. 

The side information coding portion 150 encodes the side 
information, including the information regarding the window 



employed in the adaptive block lengthMDCT 112, the scale factor 
of frequency spectrum, and the information regarding a Huffman 
table used for decoding. 

The stream composing portion 160 composes a stream of 
signals output from the Huffman coding portion 140 and the 
side information coding portion 150, with CRC (Cyclic 
Redundancy Check) added, and outputs it to the communication 
channel or the like. 

On the other hand, the decoder 200 comprises a stream 
decomposing portion 210, a scale factor decoding portion 220, 
a Huffman table decoding portion 230, a Huffman decoding 
portion 240, an inverse quantization portion 250 , and a hybrid 
filter bank 260. 

The stream decomposing portion 210 decomposes the stream 
data received from the communication channel into the frames. 

The scale factor decoding portion 220 decodes a frame 
regarding the scale factor among the data divided by the stream 
decomposing portion 210, and outputs the decoded frame to the 
inverse quantization portion 250. 

The Huffman table decoding portion 230 decodes a frame 
regarding the Huffman table among the data dividedby the stream 
decomposing portion 210, and outputs the decoded frame to the 
Huffman decoding portion 240. 

The Huffman decoding portion 240 makes the Huffman coding 
of a frame regarding the frequency spectrum of audio signals 
among the frames decomposedby the receptionprocessingportion 
210 on the basis of the information regarding the Huffman table 
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input from the Huffman table decoding portion 230, and outputs 
the decoded data (Huffman decoded data) to the inverse 
quantization portion 250. 

The inverse quantization portion 250 makes the inverse 
5 quantization of the Huffman decoded data on the basis of the 
scale factor input from the scale factor decoding portion 220 
and the Huffman decoded data input from the Huffman decoding 
portion 240, and outputs the dequantized data to the hybrid 
filter bank (HFB) 260. 
10 TheHFB 260 further comprises an alias reductionbutterf ly 

2 61, an inverse MDCT 2 62, and a sub-band synthesis filter bank 
263. 

The alias reduction butterfly 261 makes a predetermined 
butterfly operation on the dequantized data input from the 
15 inverse quantization portion 250, and outputs the operation 
result to the inverse MDCT 262. 

The inverse MDCT 2 62 that comprises an inverse MDCT buffer 
262a (see FIGS. 12 and 14 as hereinafter described) for use 
in theprocessingmakes the inverseMDCT of the operation result 
20 input from the alias reduction butterfly 261 in accordance 
with a predetermined transformation expression, and outputs 
the processed result to the sub-band synthesis filter bank. 
The inverse MDCT 262 makes a processing corresponding to each 
of the windows, including a long window and a short window 
25 or the like, on the basis of the information regarding the 
window contained in the side information. 
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The sub-band synthesis filter bank 2 63 that comprises 
a sub-band synthesis buffer 263a makes the sub-band synthesis 
of the processed result input from the inverse MDCT 262 to 
be output as an audio signal in the time domain. 
5 Now, the processing in the sub-band synthesis filter bank 

263 and the inverse MDCT 262 will be set forth. 

The sub-band synthesis filter bank 263 performs an 
arithmetical operation for an input data stream S in accordance 
with the following expression: 

10 V(i) = J^S(k)cos[x/64{i + l6X2k + l)] (o</<63) (1) 

Expression (1) canbe represented in terms of a determinant 
in the following manner. 

V = P 6 4x3 2 S (2) 

The sub-band synthesis filter bank 263 calculates an 
15 output data stream V (PCM sample data) composed of 32 data 
blocks from the input data stream S in accordance with the 
expression (2) . At this time, the sub-band synthesis filter 
bank 263 performs the multiplication by 64 x 32 = 2048 times 
and the addition by 64 x 31 = 1984 times. 
20 The input data stream S as above mentioned is the data 

rearrangedat an output signal of the inverseMDCT 2 62 . Namely, 
one time sample is extracted from each of 32 data blocks 
contained in the output signal of the inverse MDCT 2 62, and 
the extracted data of 32 samples (sample data in the frequency 
25 domain) are rearranged to generate the data block in a unit 
of 32 PCM (Pulse Coded Modulation) samples (sample data in 
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the time domain) . This processing is made for all the output 
signals from the inverse MDCT 2 62, so that the input data stream 
S composed of 18 data blocks can be obtained. 

Referring now to FIG. 11, a processing in the sub-band 
5 synthesis filter bank 263 (i.e., a sub-band synthesis 

processing) employing the above transformation expression 
will be described below. 

FIG. 11 is a diagram showing a processing procedure in 
the case where the sub-band synthesis filter bank 263 processes 
10 the input data stream S. 

In FIG. 11, the sub-band synthesis buffer 263a stores 
16 data blocks SO-0 to S0-15 calculated from the input data 
stream S at the former stage in accordance with the expression 
(2), each data block being composed of 64 pieces of data. 
15 At first, the sub-band synthesis filter bank 263 shifts 

the data blocks SO-0 to SO-15 to the right by 64 (one data 
block) (step SI) . Then, the sub-band synthesis filter bank 
2 63 makes the transformation process for the input data stream 
S in accordance with the expression (2) (step S2) , and stores 
20 a calculated data block Sl-0 in the sub-band synthesis buffer 
263a (step S3) . 

Then, the sub-band synthesis filter bank 263 extracts 
one-half (512) pieces of data from the data blocks Sl-0 to 
Sl-15 stored in the sub-band synthesis buffer 263a (step S4) . 
25 And the sub-band synthesis filter bank 263 multiplies the 
extracted 512 pieces of data by the window (step S5) , and 
periodically adds the calculated 512 pieces of data at a period 
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of 32 to output the data block (PCM sample data) composed of 
32 pieces of data (step S6) . 

Also, the inverse MDCT 2 62 performs the arithmetical 
operation on the input data stream X in accordance with the 
following expression. 

x(k) = N ]f 1 cos[- (k + 1 / 2 + N / 4Xj + 1 / 2)1 (0 < k < N - l) 

(3) 

where N is equal to 3 6 (long window) or 12 ( short window) . 
Expression (3) canbe represented in terms of a determinant 
in the following manner. 

x = P N xn/2 X (4) 

where [P NxN/2 l^ = cos^(k + V 2 + N/ 4\j + l 1 2)j 

The inverse MDCT 2 62 calculates an output data stream 
x composed of 32 data blocks from the input data stream X in 
accordance with the expression (4) . At this time, the inverse 
MDCT 262 performs the multiplication by 36 x 18 = 648 times 
and the addition by 36 x 17 = 612 times, when performing the 
arithmetical operation employing the window of the longwindow, 
or performs the multiplication by 12 x 6 = 72 times and the 
addition by 12 x 5 = 60 times, when performing the arithmetical 
operation employing the window of the short window. 

Also, theinverseMDCT2 62makesthedataprocessing (i.e., 
inverse MDCT processing) as shown in FIGS. 12 and 14 in 
accordance with the kind of window employed for processing. 



FIG. 12 is a diagram showing a processing procedure in 
the case where the inverse MDCT 262 processes the data block 
having passed through the long window. 

In FIG. 12, the input data stream X is composed of 32 
5 data blocks Dl-0 to Dl-31, each data block being composed of 
18pieces of data. Also, inFIG. 12, the datablock Dl-n denotes 
the "n" th data block among the data blocks Dl-0 to Dl-31. 
The same notation is appropriately employed below. 

At first, the input data stream X is processed through 
10 the inverse MDCT for each data block in accordance with the 
expression (3) , and transformed into the data blocks D2-0 to 
D2-31 composed of 36 pieces of data. In FIG. 12, the arrow 
indicated within the data block D2-n indicate the data 
arranging order. The same notation is appropriately employed 
15 below. 

Then, the data block processed through the inverse MDCT 
is multiplied by the predetermined window, and transformed 
into the data blocks D3-0 to D3-31 composed of 36 pieces of 
data . 

20 FIG. 13 is a diagram showing the shape of the window to 

be multiplied in the processing of FIG. 12. As shown in FIG. 
13, the windows (normal window, start window, stop window) 
may be appropriately employed in accordance with the data 
blocks D2-0 to D2-31. 

25 Of the data blocks D3-0 to D3-31, the latter half 18 pieces 

of data contained in each data block are stored in the inverse 
MDCT buffer, while the former half 18 pieces of data are added 
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to the latter half 18 pieces of data in the data block at the 
former stage, and transformed into the data blocks D4-0 to 
D4-31 composed of 18 pieces of data. 

FIG. 14 is a diagram showing a processing procedure in 
5 the case where the inverse MDCT 262 processes the data block 
having passed through the short window. 

In FIG. 14, the input data stream X is composed of 32 
data blocks Dl-0 to Dl-31, each data block being composed of 
18 pieces of data, in the same manner as in the long window. 

10 However, in the case of the short window, each of the data 
blocks Dl-0 to Dl-31 is delimited into three sets at every 
six pieces of data processed in the same short window, when 
encoded, each set of delimited data being subjected to the 
inverse MDCT processing, and multiplied by the window. 

15 FIG. 15 is a diagram showing the shape of the window to 

be multiplied in the processing of FIG. 14. As shown in FIG. 
15, the window multiplied in the processing of FIG. 14 is the 
short window, and has the shape as shown in FIG. 15. 

And each of the data blocks D5-0 to D5-31 composed of 

20 three sets of delimited data calculated in this manner is 
overlapped with and added to the data block at the former stage, 
as shown in FIG. 14, and transformed into the data blocks D6-0 
toD6-31. Moreover, the latter half 12 pieces of data contained 
in the data blocks D6-0 to D6-31 are stored in the inverse 

25 MDCT buffer, and employed for processing the data block at 
the latter stage. 
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SUMMARY OF THE INVENTION 
However, in the conventional decoder, the percentage of 
the computational complexity of the hybrid filter bank to the 
computational complexity of the entire decoder reaches about 
5 7 0 percent. Namely, in the decoder, the percentage of the 
computational complexity is greater in the sub-band synthesis 
processing and the inverse MDCT processing or the like as shown 
in FIGS. 11, 12 and 14. In such processing, the sub-band 
synthesis buffer 263a needs the amount of data of 64 x 16 = 
10 1024 blocks, and the inverse MDCT buffer needs the amount of 
data of 18 x 32 = 576 blocks, causing the cost of the entire 
decoder to be increased. 

It is an object of the present invention to provide a 
filtering method and apparatus in which the number of 
15 arithmetical operations required in the filtering process in 
decoding the digital signal and the cost of decoding is reduced . 

In order to accomplish the above object, according to 
the present invention as set forth in claim 1 , there is provided 
a filtering method for use in decoding a digital signal from 
20 the frequency domain to the time domain, including: a first 
step of multiplying an input data stream and a transformation 
matrix that is decomposed into a sparse matrix from an inverse 
MDCT transformationmatrix (e.g., a transf ormationmatrix P NxN/2 
in the expression (4)) for making the inverse MDCT 
25 transformation of the input data stream composed of a plurality 
of data blocks, and has a smaller size than the inverse MDCT 
transformation matrix, to acquire an output data stream 
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composed of a plurality of data blocks; a second step of storing 
predetermined data (e.g., a "former stage block data" as 
described in Detailed Description of the Invention) contained 
in each data block of the output data stream; and a third step 
5 of generating the digital signal in the time domain on the 
basis of each data block acquired at the first step and the 
predetermined data stored at the second step of processing 
the former stage data block. 

According to the invention as set forth in claim 9, there 

10 is provided a filtering apparatus for decoding a digital signal 
from the frequency domain to the time domain, comprising 
transformation means for multiplying an input data stream and 
a transf ormation matrix that is decomposed into a sparsematrix 
from an inverse MDCT transformation matrix for making the 

15 inverse MDCT of the input data stream composed of a plurality 
of data blocks and has a smaller size than the inverse MDCT 
transformation matrix, to acquire an output data stream 
composedof aplurality of datablocks, memorymeans for storing 
predetermined data contained in each data block of the output 

20 data stream (e.g. , an inverse MDCT buff er 12a as shown in FIGS . 
4 and 6) , and digital signal output means for outputting the 
digital signal in the time domain on the basis of each data 
block contained in the output data stream and the data of data 
block at the former stage stored in the memory means. 

25 Additionally, the invention according to claim 2 is a 

filtering method as set forth in claim 1, wherein the first 
step further comprises a fourth step of making the DCT (Discrete 
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Cosine Transform) -IV transformation process (e.g., 
transformation process in accordance with the expression (5) ) 
for each data block of the input data stream, the second step 
further comprises a fifth step of storing apart of theprocessed 
5 result of each data block at the fourth step (e.g., a "former 
stage block data" as described in Detailed Description of the 
Invention) that is used in processing the data block at the 
latter stage, and the third step further comprises a sixth 
step of folding back and expanding a part of the processed 

10 result of each data block at the fourth step other than stored 
at the fifth step and multiplying the expanded data by a 
predetermined window to acquire a first processed result, and 
folding back and expanding a part of the data block at the 
former stage that is stored at the fifth step, and multiplying 

15 the expanded data by a predetermined window to acquire a second 
processed result, and a seventh step of adding the first and 
second processed results acquired at the sixth step. 

Further, the invention according to claim 1 0 is a filtering 
apparatus as set forth in claim 9, wherein the transformation 

20 means makes the DCT (Discrete Cosine Transform) -IV 

transformation process for each data block of the input data 
stream, the memory means stores a part of the processed result 
of each data block through the DCT-IV transformation process 
that is used in processing the data block at the latter stage, 

25 and the digital signal output means further comprises: 

multiplication means for folding back and expanding a part 
of the processed result of each data block processed through 
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the DCT-IV transformation process other than stored in the 
memory means, and multiplying the expanded data by a 
predetermined window to output a first processed result, and 
folding back and expanding a part of the data block at the 
former stage that is stored in the memory means , andmultiplying 
the expanded data by a predetermined window to output a second 
processed result, and addition means for adding the first and 
second processed results output by the multiplication means . 

Further, the invention according to claim 3 is a filtering 
method as set forth in claim 1 or claim 2, wherein at the third 
step, the processing method for generating the digital signal 
in the time domain is switched on the basis of the side 
information including the information regarding the decoding 
of the data block. 

Furthermore, the invention according to claim 11 is a 
filtering apparatus as set forth in claim 9 or claim 10, wherein 
the digital signal output means switches the processing method 
for generating the digital signal in the time domain on the 
basis of the side information including the information 
regarding the decoding of the data block. 

Further, the invention according to claim 4 is a filtering 
method as set forth in claim 2, wherein at the sixth step, 
the window by which each data block is multiplied is switched 
on the basis of the side information including the information 
regarding the decoding of the data block. 

Further, the invention according to claim 12 is a filtering 
apparatus as set forth in claim 10, wherein the multiplication 
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means switches the windowby which each data block is multiplied 
on the basis of the side information including the information 
regarding decoding the data block- 
Further, the invention according to claim 5 is a filtering 
5 method as set forth in any one of claims 1 to 4, wherein at 
the first step, the output data stream is acquired by making 
the matrixing on each data block of the input data stream in 
accordance with the following expression. 

y = c-x 

10 where x is a column vector of one data block contained 

in the input data stream, y is a column vector of output data 
block corresponding to x, and C IV N is a DCT-IV transformation 
matrix represented by the following expression, 



Further, the invention according to claim 13 is a filtering 
apparatus as set forth in any one of claims 9 to 12, wherein 
the transformation means acquires the output data stream by 
making the matrixing on each data block of the input data stream 
20 in accordance with the following expression, 

y = Cx 

where x is a column vector of one data block contained 
in the input data stream, y is a column vector of output data 
block corresponding to x, and C IV N is a DCT-IV transformation 
25 matrix represented by the following expression, 




0 < m, n < N - 1 



15 



where N is the number of data contained in x. 
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[cL = 

where N is the number of data contained in x. 
Further, the invention according to claim 6 is a filtering 
method for use in decoding a digital signal from the frequency 
5 domain to the time domain, comprising a first step of 

multiplying an input data stream and a transformation matrix 
that is decomposed into a sparse matrix from a sub-band 
synthesis transformation matrix for making the sub-band 
synthesis of the input data stream composed of a plurality 

10 of datablocks andhas a smaller size than the sub-band synthesis 
transformation matrix, to acquire an output data stream 
composed of a plurality of data blocks (e.g., step S200 in 
FIG. 8), a second step of shifting each data block of the output 
data stream calculated and storedbef orehand by one datablock, 

15 every time each data block of the output data stream is 

calculated (e.g. , step S100 in FIG. 8) , a third step of storing 
the calculated data block (e.g., step S300 in FIG. 8), and 
a fourth step of generating the digital signal in the time 
domain on the basis of the stored data (e.g., step S400 to 

20 step S600 in FIG. 8) . 

The invention according to claim 14 is a filtering 
apparatus for decoding a digital signal from the frequency 
domain to the time domain, comprising transformation means 
for multiplying an input data stream and a transformation 

2 5 matrix that is decomposed into a sparse matrix from a sub-band 



f m + -Yn + -V 
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synthesis transformation matrix for making the sub-band 
synthesis of the input data stream composed of a plurality 
of datablocks andhas a smaller size than the sub-band synthesis 
transformation matrix, to acquire an output data stream 
5 composedof aplurality of datablocks, memorymeans for storing 
the calculated datablocks (e.g., the sub-band synthesis buff er 
13a in FIG. 8) , memory control means for shifting each data 
block stored in the memory means by one data block, every time 
each data block of the output data stream is calculated and 
10 storing the calculated data block in the memory means, and 
digital signal output means for outputting the digital signal 
in the time domain on the basis of the data stored in the memory 
means . 

The invention according to claim 7 is a filtering method 
15 as set forth in claim 6, wherein the first step further comprises 
a fifth step of making the DCT (Discrete Cosine Transform) -II 
transformation process for each data block of the input data 
stream, the fourth step further comprises a sixth step of 
folding back and expanding the predetermined data of the stored 
20 data block andmultiplying the expanded data by a predetermined 
window to acquire a data stream composed of a plurality of 
data blocks (e.g., step S400 and step S500 in FIG. 8), and 
a seventh step of adding the data contained in the data stream 
calculated at the sixth step at a predetermined period (e.g., 
25 step S600 in FIG. 8) . 

Further, the invention according to claim 15 is a filtering 
apparatus as set forth in claim 14, wherein the transformation 
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means makes the DCT (Discrete Cosine Transform) -II 
transformation process for each data block of the input data 
stream, and the digital signal output means further comprises 
multiplication means for folding back and expanding a 
5 predetermined data of the data block stored in the memory means , 
and multiplying the expanded data by a predetermined window, 
to acquire a data stream composed of a plurality of data blocks, 
and addition means for adding the data contained in the data 
stream that are calculated by the multiplication means at a 

10 predetermined period. 

Further, the invention according to claim 8 is a filtering 
method as set forth in claim 6 or claim 7, wherein at the first 
step, the output data stream is acquiredby making the matrixing 
on each data block of the input data stream in accordance with 

15 the following expression, 
Y = C"x 

where x is a column vector of one data block contained 
in the input data stream, y is a column vector of output data 
block corresponding to x, and C 11 ^ is a DCT-II transformation 
20 matrix represented by the following expression, 

where N is the number of data contained in x. 
Furthermore, the invention according to claim 16 is a 
filtering apparatus as set forth in claim 14 or claim 15, wherein 
25 the transf ormationmeans acquires the output data streammaking 
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the matrixing on each data block of the input data stream in 
accordance with the following expression, 
Y = C^x 

where x is a column vector of one data block contained 
5 in the input data stream, y is a column vector of output data 
block corresponding to x, and C 11 ^ is a DCT-II transformation 
matrix represented by the following expression, 

[ c «l_,„ = cos[ m ( 2n 2 + 0 < m, n < N - 1 

where N is the number of data contained in x. 

10 According to the invention as set forth in claims 1 to 

16, in decoding the digital signal from the frequency domain 
to the time domain, the number of arithmetical operations 
performed in the filtering and the memory capacity of data 
required for the operation can be reduced. Therefore, the 

15 decoding can be performed more rapidly, and the cost required 
for decoding reduced. Further, the power consumption 
required for the decoding process can be saved. 

Also, according to the invention as set forth in claims 

3 and 7 in the decoding process, a suitable decoding method 
20 can be chosen in accordance with each input data stream. 

Also according to the invention as set forth in claims 

4 and 16, with this invention, even when the coded data through 
various windows is mingled in the input data stream, the 
decoding process can be appropriately performed. 

25 

BRIEF DESCRIPTION OF THE DRAWINGS 
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FIG. 1 is a diagram showing the configuration of a hybrid 
filter bank 10 according to the present invention; 

FIG. 2 is a diagram showing the process of a matrixing 
in expression (5) ; 
5 FIG. 3 is a diagram showing the process of a matrixing 

in expression (6); 

FIG. 4 is a diagram showing a processing procedure in 
the case where an inverse MDCT 12 processes the data block 
processed in a short window; 
10 FIG. 5 is a diagram showing the shape of a windowmultiplied 

in the processing of FIG. 4; 

FIG. 6 is a diagram showing a processing procedure in 
the case where the inverse MDCT 12 processes the data block 
processed in a long window; 
15 FIG. 7 is a diagram showing the shape of a windowmultiplied 

in the processing of FIG. 6; 

FIG. 8 is a diagram showing a processing procedure in 
the case where a sub-band synthesis filter bank 13 makes a 
sub-band synthesis processing; 
20 FIG. 9 is a block diagram showing the functional 

configuration of an encoder 100 and a decoder 200 in accordance 
with the MP3 standard; 

FIG. 10 is a diagram showing a variety of windows for 
use in the coding and decoding; 
25 FIG. 11 is a diagram showing a processing procedure in 

the case where a sub-band synthesis filter bank 263 processes 
an input data stream S; 
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FIG. 12 is a diagram showing a processing procedure in 
the case where an inverse MDCT 262 processes the data block 
processed in the long window; 

FIG. 13 is a diagram showing the shape of a window 
5 multiplied in the processing of FIG. 12; 

FIG. 14 is a diagram showing a processing procedure in 
the case where the inverse MDCT 262 processes the data block 
processed in the short window; and 

FIG. 15 is a diagram showing the shape of a window 
10 multiplied in the processing of FIG. 14. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The preferred embodiments of a hybrid filter bank 10 
according to the present invention will be described below 
15 with reference to the drawings. 

FIGS. 1 to 8 are diagrams showing the hybrid filter bank 
(hereinafter referred to as an "HFB" ) 10 to which the invention 
is applied. 

The hybrid filter bank 10 according to the invention is 
20 provided in a digital signal decoder such as an MP 3 decoder, 
to transform a signal from the freguency domain to the time 
domain. 

Firstly, the configuration of the hybrid filter bank 10 
will be set forth. 
25 FIG. 1 is a diagram showing the configuration of the hybrid 

filter bank 10 to which the invention is applied. 
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In FIG. 1, the HFB 10 comprises an alias reduction 
butterfly 11, an inverse MDCT 12, and a sub-band synthesis 
filter bank 13. The configuration of the alias reduction 
butterfly 11 is the same as that of the conventional HFB 2 60, 
5 and the explanation thereof is omitted here. 

The inverse MDCT 12 comprises an inverse MDCT buffer 12a 
used in the processing (see FIG. 4 or 6 as will be described 
later) . And the inverse MDCT 12 makes an inverse MDCT in 
accordance with a predetermined transformation expression, 
10 multiplication of a window, and an overlap addition or the 
like for an operation result (a signal in the frequency domain 
that is MDCT transformed) input from the alias reduction 
butterfly 11, and outputs a processed result to the sub-band 
synthesis filter bank 13. The inverse MDCT 12 makes a 
15 processing corresponding to the window for a long window or 
short window on the basis of the information regarding the 
window contained in the side information. 

Herein, the processing in the inverse MDCT 12 will be 
now set forth. 

20 The inverse MDCT 12 multiplies an input data stream X 

that is the operation result of the alias reduction butterfly 
11 by a transformation matrix as follows. 




(5) 
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Expression (5) is derived by decomposing the inverse MDCT 
transformation matrix P n *n/2 in the expression (4) into sparse 
matrix. In the expression (5) , I N/4 is a unit matrix of N/4 
x N/4, J N/4 is an inversion matrix of N/4 x N/4, C IV is a DCT-IV 
5 transformation matrix, and "0" is a zero matrix. 

FIG . 2 is a diagram showing an operation process of matrix 
calculation in the expression (5) . As shown in FIG. 2, the 
first and second matrices on the right hand in the expression 
(5) consist of unit matrixes, inversion matrixes, and zero 

10 matrixes, and simply involve replacing and distributing the 
data. Namely, in the expression (5), the processing with the 
DCT-IV transformation matrix C IV only involves multiplication 
and addition in practice. 

Accordingly, in the case where the inverse MDCT 

15 transformation is performed in accordance with the expression 
(5) , the number of operations is halved as compared with the 
conventional method of making the inverse MDCT transformation 
in accordance with the expression (4) . For example, in a case 
where the calculation process is performed employing a window 

20 of long window, the multiplication is performed by 18 x 18 
= 324 times, and the addition is performed by 18 x 17 = 306 
times. In a case where the calculation process is performed 
employing a window of short window, the multiplication is 
performed by 6 x 6 = 36 times, and the addition is performed 

25 by 6 x 5 = 30 times. 

The inverse MDCT 12 makes the inverse MDCT for 36 points, 
thereby generating 3 6 pieces of data from the data block 
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composed of 18 pieces of data obtained by the DCT-IV 
transformation. At this time, the latter half 9 pieces of 
data among 18 pieces of data obtained by the DCT-IV 
transformation are foldedback to have 18 pieces of data, which 
5 are then directly employed for the processing with the data 
block at the former stage. The former half 9 pieces of data 
are foldedback to have 18 pieces of data, which are then employed 

y for the processing with the data block at the latter stage. 

ijj Accordingly, in processing the data block at the latter stage , 

}p 10 it is reguired to hold 18 pieces of data generated by folding 

if; 

back the former half 9 pieces of data. Herein, instead of 
storing the 18 pieces of data, the former half 9 pieces of 
data in the data block composed of 18 pieces of data obtained 
through the DCT-IV transformation are stored in the inverse 
Y 15 MDCT buffer 12a. Therefore, the amount of data stored up to 
processing the data block at the latter stage is reduced by 
half, and the memory capacity is decreased. 

The sub-band synthesis filter bank 13 comprises a sub-band 
synthesis buffer 13a. The sub-band synthesis buff er 13amakes 

20 the sub-band synthesis processing for the processed result 
input from the inverse MDCT 12, and conseguently outputs a 
digital signal (output data stream V) in the time domain. 

Herein, the processing in the sub-band synthesis filter 
bank 13 will be set forth. 

25 The sub-band synthesis filter bank 13 multiplies an input 

data stream S that is the processed result of the inverse MDCT 
12 by a transformation matrix P 6 4x32 as follows. 
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P 64x32 — L 32 C 32 (^) 

l„ 0 _ 

WherS L " =[o -co{il(2k + l)] 

Expression (6) is derived by decomposing the 
transformation matrix P 6 4x32 in the expression (2) into sparse 
5 matrix. In the expression (6), Ii 6 is a unit matrix of 16 x 
16, J 31 is an inversion matrix of 31 x 31, C 11 is a DCT-II 
transformation matrix, and "0" is a zero matrix. 

FIG . 3 is a diagram showing an operation process of matrix 
calculation in the expression (6) . As shown in FIG. 3, the 
10 first matrix on the right hand in the expression (6) consist 
of unit matrixes, an inversion matrix, and zero matrixes, and 
simply involve replacing and distributing the data. Namely, 
in the expression (6), the processing with the DCT-II 
transformation matrix C 11 only involves multiplication and 
15 addition in practice. 

Accordingly, in the case where the sub-band synthesis 
processing is performed in accordance with the expression ( 6) , 
the multiplication is performed by 32 x 32 = 1024 times, and 
the addition is performed by 32 x 31 = 992 times, whereby the 
20 number of operations is halved as compared with the 

conventional method making the sub-band synthesis processing 
in accordance with the expression (1) . 

Moreover, in the sub-band synthesis processing, it is 
required to hold 64 pieces of data that is the multiplication 
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result of the input data stream by the transformation matrix 
in the expression (6) to be used in the data block at the latter 
stage. However, these 64 pieces of data are obtained by 
replacing and distributing 32 pieces of data obtained by the 
5 DCT-II transformation. Accordingly, instead of storing the 
64 pieces of data, the 32 pieces of data are stored in the 
sub-band synthesis buffer 13a, whereby the memory capacity 
is reduced by half as compared with the conventional method. 
The operation will be described below. 
10 First of all, the operation of the inverse MDCT 12 will 

be set forth. 

FIG. 4 is a diagram showing a processing procedure for 
the inverse MDCT 12 to process the data block that has passed 
through the short window. 

15 In FIG. 4, an input data stream is composed of 32 data 

blocks Dl-0 to Dl-31, each data block being composed of 18 
pieces of data. Each data block is delimited into three sets 
for every 6 pieces of data having passed through the same short 
window, when encoding, and the inverse MDCT transformation 

20 (in accordance with the expression (5) ) is made for each 
delimited data (hereinafter referred to as a "short block") . 
In FIG. 4, the notation "x (-1) " means that each value of data 
array is multiplied by "-1". The same notation is 
appropriately employed in the following. 

25 Three short blocks having undergone the inverse MDCT 

transformation processing (hereinafter referred to as the 
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first to third short blocks in the order of short block array, 
each designated with © to (3) in FIG. 4) . 

Firstly, the entire data in the first and second short 
blocks and the latter half of data in the third short block 
5 (latter three pieces of data) are folded back and expanded 
into double size. 

Also, the 9 pieces of data in the data block at the former 
stage (hereinafter referred to as a "former stage block data") 
that are stored in the inverse MDCT buffer 12a, and expanded 

10 into four-thirds times (12 pieces of data) . These pieces of 
data are multiplied by respective windows, and the 
multiplication results are over lapping addition, as shown 
in FIG. 4. Consequently, the data block L composed of 24 pieces 
of data is generated. 

15 Herein, the former stage block data is processed by 

switching the window depending on whether the data is encoded 
through the long window or the short window and so on. 

FIG. 5 is a diagram showing the shape of window that is 
multiplied in the processing of FIG. 4. As shown in FIG. 5, 

20 in the case where the former stage block data is encoded through 
the long window, the former stage block data is multiplied 
by the latter half 18 pieces of data in the start window. Since 
the last 6 pieces of data in the start window are "0", the 
description is omitted in FIG. 5. Also, in the case where 

25 the former stage block data is encoded through the short window, 
the former stage block data is multiplied by the latter half 
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6 pieces of data in the short window. Moreover, each short 
block is multiplied by the short window in either case. 

Among the data contained in the generated data block L, 
the former half 18 pieces of data are output as the operation 
5 result of the inverse MDCT 12. On the other hand, the latter 
half 6 pieces of data contained in the data block L are connected 
with the latter half 3 pieces of data in the third short block, 
as above described, and stored at a predetermined address of 
the inverse MDCT buffer 12a. 

1° FIG. 6 is a diagram showing a processing procedure for 

the inverse MDCT 12 to process the data block that has been 
processed by the long window. 

In FIG. 6, an input data stream is composed of 32 data 
blocks Dl-0 to Dl-31, each data block being composed of 18 

15 pieces of data. 

And the data blocks are subjected to the DCT-IV 
transformation processing (in accordance with the expression 
(5)), and transformed into data blocks D7-0 to D7-31. 

Next, each of the data blocks D7-0 to D7-31 subjected 

20 to the DCT-IV transformation processing is divided into the 
former half 9 pieces of data and the latter half 9 pieces of 
data in each data block. The latter half 9 pieces of data 
are folded back and expanded into double size together with 
the former stage data block stored in the inverse MDCT buffer 

25 12a, respectively. These pieces of data are multiplied by 
respective windows, and the multiplication results are over 
lapping addition, and output as the operation result of the 
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inverse MDCT 12. And the former half 9 pieces of data are 
stored at the predetermined address of the inverse MDCT 12. 

Herein, the former stage data block is processed by 
switching the window depending on whether the data-is encoded 
5 through the long window or the short window. 

FIG. 7 is a diagram showing the shape of window that is 
multiplied in the processing of FIG. 6. As shown in FIG. 7, 
in the case where the former stage block data is encoded through 
the normal window, the former stage block data is multiplied 

10 by the latter half 18 pieces of data in the normal window, 
and the data (referred to as the data blocks D8-0 to D8-31) 
expanded from the latter half 9 pieces of data in the data 
blocks D7-0 to D7-31 are multiplied by the former half 18 pieces 
of data in the normal window. Also, in the case where the 

15 former stage block data is encoded through the short window, 
the former stage block data is multiplied by the latter half 
6 pieces of data in the start window, and the data blocks D8-0 
to D8-31 are multiplied by the former half 18 pieces of data 
in the stop window. 

20 As shown in FIGS. 4 and 6, the processed result by the 

inverse MDCT 12 in which the window is switched in accordance 
with the processing object data block is multiplied by the 
transformation matrix of the expression (6) in the sub-band 
synthesis filter bank 13, and output as a digital signal in 

25 the time domain. 

The operation of the sub-band synthesis filter bank 13 
will be now set forth. 
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FIG. 8 is a diagram showing a processing procedure for 
the sub-band synthesis filter bank 13 to make the sub-band 
synthesis processing. 

In FIG. 8, 16 data blocks S2-0 to S2-15 transformed from 
5 the input data stream S at the former stage in accordance with 
the expression (6) are stored in the sub-band synthesis buff er 
13a, each data block being composed of 32 pieces of data. 

At first, the sub-band synthesis filter bank 13 shifts 
the data blocks S2-0 to S2-15 stored in the sub-band synthesis 
10 buffer 13a to the right by one data block (step S100) . 

Also, the sub-band synthesis filter bank 13 makes the 
32 point DCT-II transformation process for the input data 
stream S for each data block to calculate the data block S3-0 
successively (step S200) . And the sub-band synthesis filter 
15 bank 13 stores one calculated data block S3-0 in the sub-band 
synthesis buffer 13a with the data shifted by one data block 
(step S300) . 

Then, the sub-band synthesis filter bank 13 rearranges 
the 32 x 16=512 pieces of data stored in the sub-band synthesis 
20 buffer 13a in a predetermined order (step S400) . 

Then, the sub-band synthesis filter bank 13 multiplies 
the 512 pieces of rearranged data by the window, and outputs 
16 data blocks composed of 32 pieces of data (step S500) . 
Lastly, the sub-band synthesis filter bank 13 adds 16 
25 data blocks periodically, and outputs the data block (PCM 
sample data) composed of 32 pieces of data (step S600) . 



T 
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The sub-band synthesis filter bank 13 outputs the PCM 
sample data by repeating step S100 to step S300 successively 
for every processing of one data block contained in the input 
data stream. 

5 As described above, with the hybrid filter bank 10 to 

which this invention is applied, in performing the inverse 
MDCT process, the inverse MDCT 12 makes the transformation 
process in accordance with the expression (5) . 

Accordingly, the number of arithmetical operations made 

10 in the inverse MDCT transformation process is reduced by half 
as compared with the conventional method, resulting in the 
faster processing. 

The former stage block data stored in the inverse MDCT 
transformation process is composed of 9 pieces of data. 

15 Accordingly, the memory capacity is reduced by half as compared 
with the conventional method in which it was composed of 18 
pieces of data. Because the memory capacity provided for the 
hybrid filter bank is halved, the costs can be reduced. Also, 
when the inverse MDCT transformation process is performed, 

2 0 the power consumption for storing the data can be saved. 

Moreover, the sub-band synthesis filter bank 13 performs 
the transformation process in accordance with the expression 
(6) . Accordingly, the number of arithmetical operations made 
in the sub-band synthesis process is reducedby half as compared 

25 with the conventional method, resulting in the faster 

processing. Moreover, the memory capacity reguired in the 
sub-band synthesis process is reduced by half as compared with 
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the conventional method because the intermediate operation 
results of sparse matrix decomposition are stored, resulting 
in the lower costs . Also, when the sub-band synthesis process 
is performed, the power consumption taken for storing the data 
5 can be saved. 

The hybrid filter band 10 to which the invention is applied 
makes the inverse MDCT transformation process by switching 
the window appropriately for each block of the input data stream 
on the basis of the side information. Accordingly, even when 
10 the data blocks encoded by various windows are mingled in the 
input data stream, the decoding process can be appropriately 
performed. 

This invention can be applied to the filter bank of the 
Daubechies Wavelet Transform that is employed in a still image 
15 compression method using a wavelet transformation. 

According to the present invention as set forth in claims 
1 to 16, in decoding the digital signal from the frequency 
domain to the time domain, the number of arithmetical 
operations performed in the filtering and the memory capacity 
20 of data required for the operation can be reduced. Therefore, 
the decoding can be performed more rapidly, and the cost 
required for decoding reduced. Further, the power 
consumption required for the decoding process can be saved. 
Also according to the invention as set forth in claims 
25 3 and 7, in the decoding process, a suitable decoding method 
can be chosen in accordance with each input data stream. 
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Also according to the invention as set forth in claims 
4 and 8, even when the coded data through various windows is 
mingled in the input data stream, the decoding process can 
be appropriately performed. 



